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Metal-insulator transitions and transitions between different quantum Hall liquids are used to 
describe the physical ideas forming the basis of quantum phase transitions and the methods of 
application of theoretical results in processing experimental data. The following two theoretical 
schemes are discussed and compared: the general theory of quantum phase transitions, which has 
been developed according to the theory of thermodynamic phase transitions and relies on the concept 
of a partition function, and a theory which is based on a scaling hypothesis and the renormalization- 
group concept borrowed from quantum electrodynamics, with the results formulated in terms of flow 
diagrams. 
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1. INTRODUCTION 

Quantum phase transitions that occur in an electron gas and result in localization represent the mainstream in 
the study of electrons in solids; it is directed toward low temperatures and interactions. We think that a dangerous 
gap between theory and experiment has appeared in this field of science. The related theories are so complex that 
experimentalists, as rule, use ready recipes and formulas of these theories without knowing (and, hence, without 
controlling) all the initial assumptions and limitations. 

In this review, we will try to close this gap. Our review does not contain a sequential description of the mathematical 
techniques involved in formulating hypotheses and theories. Instead, we outline the physical ideas, concepts, and 
assumptions that are often omitted in theoretical works and reviews, especially at the stage of a developed theory. 
This specific feature is thought to make our review interesting and useful for experimentalists. 

On the other hand, we did not tend to accumulate, classify, and estimate the huge experimental data obtained in 
this field. We address experiment to demonstrate the use of a theory for its interpretation and the related difficulties 
and problems. This specific feature is thought to make our review interesting and useful for theorists. 

We first qualitatively discuss the ides that constitute the basis of the theory of quantum phase transitions. We 
then consider how the general theoretical scheme can describe the well-known data on metal-insulator transitions and 
transitions between different quantum Hall liquids. 

1.1. Thermal and quantum fluctuations 

Before discussing quantum phase transitions in essence, we briefly recall some facts from statistical physics. We 
consider a macroscopic system consisting of a huge number of particles. This system is almost isolated and behaves 
mainly as a closed system; however, it can exchange particles and energy with a larger system, which serves as a 
reservoir. In other words, our system is a subsystem of this reservoir. 

We first analyze the classical subsystem at a finite temperature T. According to classical statistics formulas, the 
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probability pi for the subsystem to occupy the state with energy e$ is 

Pi = , (1) 

where the coefficient 1 jZ is determined by the normalization condition Y^Pi = 1 

Z = ^exp(- £i /T) (2) 

i 

(hereafter, temperature is given in energy units). Function Z is called the partition function. It plays a key role in 
the description of the thermodynamic properties of classical objects. That the subsystem energy obeys distribution 
([I]) rather than being fixed means the presence of classical thermal fluctuations. 

When quantum mechanics appeared, the basis of classical thermodynamics and its principal formulas were revised. 
The classical expressions were found to have a limited field of application. Quantum statistics requires a developed 
quantum- mechanics math-based environment. Before addressing this environment, we assume that the subsystem 
is completely closed and does not interact with the large system. As a result, we can use the standard quantum 
mechanics technique and use the Schrodinger steady-state equation to describe the subsystem. The resulting set of 
steady-state energies £j and the corresponding wavefunctions ipt (q) of the subsystem are considered to be subsystem 
attributes in the zeroth approximation. These energies, £j, enter Eqn ((2]). The set of <fi(q) functions is convenient 
because it is related to the subsystem under study and is a complete set; therefore, it can be used for expansion into 
a series. 

The presence of a huge number of particles in the subsystem implies that it has a very dense energy distribution of 
quantum levels. Due to a weak interaction with the reservoir, the subsystem is in the so-called mixed state, in which 
no measurement or a set of measurements can lead to unambiguously predicted results. Upon determining the mixed 
state, the repeated measurements of any physical quantity give a near-average result that differs from the previous 
one. The scatter of the measured physical quantities is interpreted as a result of fluctuations. 

Because the subsystem interacts with the environment, it cannot be described by a wave function; therefore, the 
wave function is substituted by a density matrix, 

pW,q)= I r(q',X)i;(q,X)dX, (3) 



where q corresponds to the set of subsystem coordinates, X are the remaining coordinates of the system, and ^(q, X ) 
is the wave function of the closed system. 

The mean values of any physical quantity (s) are now calculated using the density matrix 

{S) = \\ {' [P(9 '' q)] } g '= q dq ' A = J P{q ' q)dq - (4) 

rather than a wave function. In Eqn Q, we apply the operator s, which acts on functions of the variable q, to the 
density matrix p(q',q); we then set q' = q and integrate. The operator that formally enters the expression for the 
normalization factor A is identical with unity; therefore, the condition q ' = q is naturally taken into account in the 
integrand. 

According to general ruleQ, the mean of coordinate (q), for example, is given by 

M-^/^™*^ («) 

The physical meaning of the density matrix can be clarified if we write it in an explicit form for the subsystem in 
statistical equilibrium at a finite temperature T; 

p(q' } q)=Y / tf(q')^(q)e^ /T . (6) 

i 

For the subsystem in statistical equilibrium, it follows from Eqns © and ((U 

^ = J/Zj tf(qMq)e- E ' /T qdq = i^^e-^, (t? ). = J tf( q )^( q ) qdq . (7) 
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The averaging with the density matrix performed in Eqn |(7|) account for both the probabilistic description in the form 
of (q)i in quantum mechanics and incomplete information about the system (statistical averaging). 
We use the (fi(q) functions and write the function p(q', q) in the matrix form 



p(q',q) = \\pij\\ 

In Eqn (j4]), we then have 



i pUq)p(q',q)<fi(q)dq 



(8) 



/ v 12, ; S ijPji \ - Pij ,_, 

(s) = = SijWjh Wij = (9) 

Z^ii P™ , 2-ii Pit 

The Wij matrix constructed from the set of <fii{q) functions is called the statistical matrix and, in essence, is the 
normalized density matrix. 

We replace s with the energy operator H in Eqn and obtain 

(s) = ^Wu£i. (10) 

i 

This means that the probability pi of detecting the energy e, in the subsystem is equal to the diagonal element of the 
matrix wu, 

Pi = wu. (11) 

Equation l[TTj) is the quantum analog of Eqn (U) . 

The statistical matrix has a number of universal properties. By definition, this matrix is normalized: 

= 1. (12) 

i 

As follows from Eqn ((6|), the statistical matrix is diagonal in statistical equilibrium. Its diagonal elements are functions 
of only the energy of the corresponding subsystem state £j. Until at least one wu with i ^ has a nonzero value, the 
subsystem state remains mixed, and any measured quantity undergoes fluctuations. 

The macroscopic system is in a mixed state mainly due to a finite temperature. A comparison of Eqns {TJ and l[TTj) 
demonstrates that in the high-temperature limit, the statistical matrix obeys the Gibbs distribution 

w u = exp(-£i/T)/Z. (13) 

with a good accuracy. 

Because the integral of the squared wave function <pi (q) over the entire coordinate space is 

Pi(q)<Pi(q) d q = ^ 

the partition function Z can be represented as the sum of the diagonal matrix elements of the operator exp (-H/T) 
over the complete set of eigenfunctions ipc 

Z = E / V*i(q)vMe- ei/T dq = £>i| exp(-H/T)\ Vi ). (14) 

i J i 

Here, we use that ipi are eigenfunctions of the Hamiltonian H to obtain 

exp(-H/T)<pi = e- ei > T <pi 



A comparison of Eqns {!])— with Eqn lfT3|) gives another form of the partition function, 
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Z = I p(q,q)dq = S V [p(q',q)]. (15) 

If Z is regarded as a normalization factor, Eqns |(4]) and iJH) demonstrate that representation lfl~5|) is valid in both the 
classical limit and general case. 

As the subsystem is open, it interacts with the environment. Even if this interaction is very weak, the states of the 
macroscopic system become mixed in a certain small energy range due to the extremely high density of the system 
energy levels. Therefore, even at zero temperature, the unclosed subsystem is in a mixed state. This means that as the 
temperature decreases, each wu (i ^ 0) element approaches a temperature-independent constant, which depends on 
Ej and the sizes and properties of the subsystem, rather than tending to zero exponentially. Thermal fluctuations are 
replaced by quantum fluctuations. This statement is illustrated by Fig. UK, which qualitatively shows the temperature 
dependence of the probability pt of detecting the energy £j in the subsystem. The dashed lines separate the regions of 
predominantly quantum and predominantly thermal fluctuations. The limiting value Pi(0) is specified by the energy 
Ei, sizes, and properties of the subsystem. 




FIG. 1: (a) Low-temperature behavior of the diagonal elements wu of the statistical matrix corresponding to low energies, (b) 
Simultaneous excitation of thermal and quantum fluctuations 



In this subsystem, pi(0) obviously decreases with an increase in £,-, We assume that this dependence is a power-law 
function, pi(0) oc l/e", and obtain the lower estimate of the exponent a. Energy for thermal fluctuations is taken 
from the reservoir. A finite, not exponentially small, probability of detecting the subsystem in a state ei~^> T means 
a violation of the energy conservation law. This violation can exist in the framework determined by the uncertainty 
relation 



EiTi ~ h, (16) 

where Tj is the lifetime of the state &j. Therefore, even if transitions to all e, states were equiprobable, Pi(0) would 
decrease according to the law Pi(0) oc e" 1 , as follows from Eqn lfT6|) . As a result, we have 

a Ss 1. (17) 

It is essential that the energy range in which quantum fluctuations occur is temperature independent, which follows 
from uncertainty relation l|16p . 

The relative role of thermal and quantum fluctuations is illustrated in Fig.QJ), which shows the excitation spectrum 
of the subsystem in one of the phases, i.e., far from phase transition point x = x c . The brace in the bottom left 
part of Fig.QJ) indicates the temperature-related scale. Only modes with frequencies u>i < T/h are classically excited 
modes. Although the range of thermal modes is bounded above by temperature, they have large occupation numbers. 
Modes with u>i > T/h are mainly excited due to quantum processes and have small occupation numbers. However, 
as the temperature decreases, the role of quantum excitations increases. 
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2. QUANTUM PHASE TRANSITIONS 
2.1 Definition of a quantum phase transition. 

Quantum phase transitions are phase transitions that can occur at the absolute zero T — and consist of a change 
in the ground state of a system in the case where a certain control parameter x takes a critical value x c . (The ground 
state is taken to be the lowest possible mixed state.) The control parameter can be, for instance, a magnetic field or 
an electron concentration. Quantum phase transitions belong to the class of continuous transitions at which none of 
the physical functions of state has a discontinuity at the transition point. 

A quantum fluctuation is the only cause that can change the ground state of the system at zero temperature; 
the phase transitions under study are therefore called quantum phase transitions. As in the case of thermodynamic 
transitions, the concept of a correlation length £, which has the meaning of the average quantum-fluctuation size, 
is introduced into the theory of quantum phase transitions. At the absolute zero temperature, £ is only determined 
by the deviation Sx of the control parameter from the critical value. The dependence of £ on Sx is assumed to be a 
power-law function, 



tK\x-x e \- v . (18) 

Real experiments are always performed at finite temperatures, where thermal fluctuations exist in addition to quantum 
fluctuations. The goal of the theory is to predict the manifestation of the phase transition that occurs only at zero 
temperature in the properties of the system at a finite temperature. 

The theory of quantum phase transitions (see, e.g., book [l| or reviews 0, 0|) is analogous to the theory of 
thermodynamic phase transitions. In the (x,T)) plane (i.e., in a plane where temperature is plotted versus a control 
parameter), the quantum transition point can be a finite point on the line of thermodynamic transitions, x c (T) — > 
x c (0) as T — > 0. For example, this behavior is characteristic of magnetic transitions with the magnetic field used as a 
control parameter. Quantum phase transitions of another type also exist; they are imaged by isolated point x c on the 
abscissa axis of the {x,T) plane. The metal-insulator transition is an example of such a transition. In this review, 
we restrict ourselves to the case of isolated-point transitions. 



2.2 Parallels and differences between classical and quantum phase transitions 

As noted above, quantum phase transitions belong to the class of continuous transitions; there are no two coexisting 
competitive phases and, hence, no stationary boundary between them. This class also includes thermodynamic second- 
order phase transitions. At a thermodynamic phase transition point, the system transforms into another phase as a 
whole as a result of thermal fluctuations. Fluctuations exist on either side of the transition, and their characteristic 
size £ is called the correlation length. As the transition is approached from either side, £ diverges 0, 0|. 

It is natural to expect a similar situation to occur in the vicinity of a quantum transition with the participation 
of quantum fluctuations. An analogy between classical and quantum phase transitions does exist, and it is rather 
unexpected. The behavior of a quantum system in the vicinity of a transition point at a finite temperature in a d- 
dimensional space is analogous to the behavior of a classical system in a space with dimension V > d. This statement 
requires extensive explanations, which should include an algorithm for the introduction of such an imaginary system 
and the determination of its dimension. 

First, we have to clearly distinguish between the dimension of the geometric space in which the system is located 
and the dimension of the generalized-coordinate space, which depends on the number of particles N in the system. If 
q is the density of particles in the space, we have J\f = f f j gd d X , where the dimension d determines the multiplicity 
of the integral. Usually, a d-dimensional space is assumed to be infinite in all directions, but the range of one or 
several Xj coordinates can be limited. 

Second, we note a similarity between the operator exp (-H/T), which was used to write Eqn (HI), and the operator 
S that describes the evolution of a closed quantum system with time according to the Schrodinger equation 

ih^=Hip, S = exp(-~Ht). (19) 

This similarity becomes obvious if we change the variables as 

it/h = 1/T. (20) 

With substitution (|20| . we can interpret the matrix elements in Eqn (fl4| differently. The element 
(i\exp(-H/T)\i) can be regarded as the amplitude of the probability that, starting from state (i |, the subsys- 
tem evolves under the action of S and returns to the initial state \i } within imaginary time t, which is equal to —ih/T. 
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The imaginary time i is often called the Matsubara time. For definiteness, we assume that the number of steps is 
fixed and equal to N + 1, N ^> 1 and that, in time ih/T, the subsystem has passed through N virtual states and 
occupied each state for a time (St)j; as a result, we have 



JV 

= ih/T. (21) 

The amplitude of 'the probability of returning to the initial state' means the sum of the amplitudes of the proba- 
bilities of returning for all possible trajectories in the space of states. We consider a set of trajectories consisting of 
N steps. The operator entering each matrix element in Eqn lfl4|) is represented as 



N 



e-xp(-H/kT) = cxp i^H^2(5t)jj . (22) 

We take a complete set of wavefunctions \rrij) of an arbitrary operator that does not commute with the Hamiltonian 
and assume that the trajectories are realized in these states. Then, instead of Eqn l[T4|) . we obtain 



Z = ^2 ^2 <i \ cx P(-i#W)iA)l TO i >< mi\exp(-iH(St) 2 /h)\m 2 > ... < m N \ exp(-iH (St) N / h)\i > . 

i mi,m2,...mw 

(23) 

The product of the matrix elements in the summand corresponds to a chain of consecutive virtual transitions. The 
summation over all {rrij} combinations means that all possible closed chains of N links are taken into account. The 
modification of the expression for Z — the replacement of Eqn 1|14[1 with Eqn lf23|) — means that we take the quantum 
properties of the system into account by adding virtual transitions to real ones. 




FIG. 2: Set of quantum statistical d-dimensional systems located in a one-dimensional band of width ft/Tand virtual transition 
between them 



The introduction of virtual transitions is illustrated in Fig. [21 where the abscissa axis stands for the initial d- 
dimensional space { d X} and all possible classical states of the system are located along this axis in the order of 
increasing energy. All other horizontal lines are replicas of this space. These replicas form a set of N elements in the 
segment [0,h/T\. The black dot on the abscissa axis represents the initial state of the system with energy The 
second black dot on the upper horizontal line corresponds to the 'final' state of the system with the same energy e$ 
at the maximum value of the imaginary time \t | = %/T. The arrows indicate the chain of virtual transitions through 
the set of virtual states \rrij) designated by white dots. These are mixed states without a definite energy; therefore, 
they are not eigenstates of the operator H. The virtual transitions reflecting the quantum properties of the statistical 
system under study are shown by arrows. The chain of virtual transitions shown in Fig. [2] corresponds to one term in 
the sum in (f23|) . The summation over all {rrij} combinations means that we account for all possible chains between 
the fixed initial and final £j points. 

Now, the problem is to construct a classical system whose partition function is represented by Eqn (f23]) . In the 
original Eqn l[T4j) . the summation over i meant the summation over all the states realized in a d-dimensional system. 
The number of terms in the sum is now increased, and our real system does not have such a large number of different 
states. Nevertheless, we can construct an imaginary classical system with a higher dimension using the scheme shown 
in Fig. H 
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We add an additional axis for an imaginary time to the d axes of the original space. In the graphical terms of Fig. 
[21 this means that an ordinate axis is added to the abscissa axis. It is seen from Eqn (plj) that the coordinate along 
this new axis changes in the range 



< \t\ < h/kT. (24) 

Each point in the strip (|24| in Fig. [2] corresponds to some virtual state of the quantum system; it can be named the 
image point. The width lj24|) of strip in which an image point is located increases with decreasing the temperature. 
At T = 0, the strip transforms into a half-plane; an increase in the temperature, in contrast, narrows the band and 
decreases its contribution to the statistical properties of the quantum system. 

We repeat the d-dimensional classical system N times by placing replicas along the imaginary time t axis at a 
distance (St)j from each other. The states of the replicas in this ensemble can be different. Let these states be rrij. 
We fix a certain 'initial' state of the 'lower' subsystem £j. The energy of the classical 2?-dimensional system, we are 
constructing, is equal to the sum of the energies of all layers, with allowance for the interaction between them. Each 
£i,{m} corresponds to a certain energy e, of the initial d-dimensional classical system and a certain chain of states 
entering Eqn (|23l> . Therefore, the number of terms in the sum J2{ m } s i-m hi Eqn ([23| is equal to the number of 
different £jx m } energies. We note that each Sj im term in sum has a form typical for partition sums, since the product 
of the matrix elements in this term eventually reduces to the product of exp (— iEj(8t )j/n). The only difference is 
that for the classical system, each term in Z is a real positive number and, in Eqn (J23J) , a real positive number is 
represented by the sum Sj = J2{ m } s i-™ over { m } °f a U s i- m products. When neglecting this difference, we can 
consider Eqn (J23j) the partition function Z of the imaginary ©-dimensional classical system. 

Of course, in the general form, the resulting mathematical construction is absolutely unpractical. However, as is 
customary in deriving scaling relations, it is important to reveal some significant properties of sum (|23|) rather than 
to calculate it. As is shown in Section 2.3, one of these properties is anisotropy, i.e., non equivalence of the axes of 
the D-dimensional space. 



Ac 




Phase 1 



Phase 2 



FIG. 3: The difference in the lowest energies of two phases Ae = eo2 — £oi versus a control parameter (dashed line). At a 
control parameter x = x c — 8x, phase 1 is subjected to fluctuation excitation inside equilibrium phase 2. Because |Ae| < T, 
thermal fluctuation play a key role in this excitation 



At a thermodynamic phase transition point, the partition function Z has special features. The sensitivity of function 
(f!4|) to the presence of a phase transition is caused by the fact that, when the transition is approached, the energy 
range T determining significant terms in sum (f!4|) contains not only ipi levels from the set corresponding to the 
equilibrium phase but also ip\ levels from the nonequilibrium phase (Fig. EJ. This allows fluctuation transitions 
between \i } and \i '} levels from different sets. Additional possibilities appear after Eqn ifFf) is replaced with Eqn ([23)) . 
If the temperature is low (|Ae| > T), the fluctuation- induced appearance of another phase is also possible, but due 
to quantum fluctuations. 



2.3. Critical region of a quantum phase transition 

A region in which all physical quantities depend only on the correlation length £ always exists in the phase plane 
near a classical phase transition. It is called the critical or scaling region. Near a quantum phase transition at 
T = 0, there also exists a control parameter Sx range in which physical quantities are expressed through the length 
£ determined by Eqn 1(18)1 . In Section 4, we show this behavior by the example of a metal-insulator transition in 
a three-dimensional system of noninteracting electrons, and determine the boundaries of this region. At a finite 
temperature T ^ 0, however, the scaling region of a quantum phase transition is more complex. 
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The space { d X, t } has dimension d + 1 because of the additional imaginary time axis. We introduce correlation 
lengths in the space { d X,t}. We retain the traditional notation £ for the correlation length in an ordinary d- 
dimensional subspace and let £ v denote the correlation length along the additional axis. The subscript <p is a reminder 
that £ v is related to the quantum aspect of the problem and to the specific features of wave functions. The dimension 
of !; v is H/T rather than length; that is, it is measured in seconds rather than centimeters. 

As x — y x c and T — > 0, both correlation lengths diverge at a transition. According to the theory of continuous 
thermodynamic transitions, the divergence is described by power functions; but the exponents of the two correlation 
lengths can be different in general. This is usually written as 

£ oc 5x~ v , 8x=\x — x c \, (25) 



^ oc C 2 . (26) 

The exponents v and z are called critical indices; z is called the dynamic critical index. Both names and designations 
originate from the theory of thermodynamic processes. In particular, the dynamic critical index in the theory of 
thermodynamic processes enters the relation between the lifetime and size of thermal fluctuations, r oc £ z . In 
the quantum problem at T = 0, the situation is similar: the length £ characterizes spatial correlations, i.e., the 
characteristic size of quantum fluctuations, and the time £ v characterizes time correlations. 

Formula (|26l> immediately fixes the dimension of the space in which the imaginary classical system should be placed. 
A real length must be the coordinate along all axes of this space; that is, length axes should be made from the t axis. 
As follows from dimensional considerations and Eqn l(26|) . the length equivalent to the 'correlation pseudolength' £ v 
is proportional to £<J/ Z , 

Q^ej z - (27) 

Therefore, the volume element in the space is 

{d£) d d£ v <x {dO d {dCf (28) 

and a ^-dimensional subspace with the conventional 'spatial' coordinates appears instead of the one-dimensional 
imaginary time axis. Hence, we have 

V = d+z. (29) 




FIG. 4: Shape of the critical region of a quantum phase transition depicted by an isolated point in the {x,T) plane 

It is convenient to discuss the consequences of Eqns (|25| and f26|) using the (x,T) plane (Fig. SJ. At T = 0, the 
control parameter x on the abscissa axis of the (x,T) plane only affects one independent correlation length, £, and 
£ v (and the length () can be formally obtained from £ using Eqn lj26j) . All the physical quantities can be expressed 
only in terms of x c in a certain segment [xi,X2] that contains the point x c . 

Now, let the temperature T ^ 0. We move in the (x, T ) plane along a horizontal line (Fig. |4j line aa) in the 
x — -> x c direction. At a certain value of 5x, the length which changes according to Eqns l(25|) and (|26"ll . reaches its 
maximum value H/T [which is determined by inequality lf24j) ]. At lower values of Sx, Eqn lj26j) does not hold on line 
aa, because £ increases according to Eqn lj25jl and £ v remains equal to its maximum value H/T. The parameters £ 
and £ v become mutually independent. The region of this independence is called the critical region. If we now move 
toward the transition inside the critical region (e.g., along line bx c ), both parameters still diverge at the transition. 
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However, the divergence of one of them is controlled by Sx, £ oc Sx ~ v , and the divergence of the other is controlled 
by temperature, £ v = h/T . 

The name of the region derives from the fact that, for a quantum transition, the critical region is considered to be 
the region where quantum fluctuations play an essential role in mixing the two-phase states. As can be seen from 
Fig. [3l this occurs only under the condition 

|e 2 -soil > T. (30) 



As the temperature decreases, the range of the control parameter Sx where condition ([30)1 is satisfied narrows. As a 
result, the shape of the critical region is unusual: it widens when moving from the transition. 

Inside the critical region, the quantity £ v = H/T can be associated with a real length. The quantity 

L v oc (h/T) 1 /*, (31) 

has the required properties: it is equivalent to the parameter £ v = h/T and, according to Eqn lj26j) . has the dimension 
of length. The physical meaning of L v can be understood from the following considerations. Because of a finite tem- 
perature, the quantum problem acquires a characteristic energy T that separates classical and quantum fluctuations 
(see Fig. [TJd). Quantum fluctuations are fluctuations with energies hu v > T. The spatial size of these fluctuations 
l v oc \/lo v is bounded above because frequency lo v > T/h is bounded below. The length L v is the upper boundary of 
the size of quantum fluctuations. Therefore, L v is often called the dephasing length, i.e., the length beginning from 
which coherence in a system of electrons is broken. 

The appearance of two independent parameters in the critical region of a quantum phase transition is caused by 
inequality ([24| . Therefore, a scaling description in the critical region of a quantum phase transition is therefore called 
a finite-size scaling. An imaginary thermodynamic system is thought to exist in a hyperstrip in a D-dimensional space 
with d variables ranging from to oo and z variables ranging from to (h/T) 1 /*. 

Inside the critical region, L v < £, and, along its boundary, 

L v = e (32) 

The lengths L v and £ are determined by Eqns ([25[) and ([26j) up to a multiplicative constant; however, Eqns ([25| and 
([26[) rigidly fix a power relation between these variables. Therefore, the equation for the critical-region boundaries 
has the form 

T = C{5x) vz , (33) 

where the constant C can be different on either side of the transition in general. 

Because the critical region has two independent parameters, the scaling expressions for physical quantities become 
more complex. We only present and discuss expressions for the electric conductivity and resistivity, 

a = e- d F(L v /0, p = Z d - 2 F 1 (L v /0, F 1 {u) = l/F{u), (34) 

where F(u) is an arbitrary function. The exponent of the first factor in Eqn lj34|) is specified by how the length enters 
the expressions for conductivity at various dimensions d. As can be seen from Fig. [4j this factor is determined by 
the z-component of the distance from the transition. The L v /£ ratio is dimensionless; therefore, scaling rules can be 
used to regard L v /£ as an argument of an arbitrary function. This relation depends on both T (i.e., the y-component 
of the distance from an image point to the transition in the phase plane) and the distance from this point to the 
critical-region boundary along the x-axis. 

As an argument, we can also take any power of the L v /^ ratio; as a result, we write the argument in various forms, 
such as 

h/T_ H/T h{5x) z » Sx 

These various forms elucidate various aspects of the physical meaning of this ratio. In the last two forms, the 
argument is dimensional and the dimensionality of the control parameter is not specified a priori. These forms show 
how temperature enters the argument of the scaling function. 

In conclusion of this section, as an example, we present an implication of Eqn ([34| for the resistivity of a two- 
dimensional system. If a quantum phase transition induced by electron localization occurs in a two-dimensional 
system and manifests itself in the transport properties of this system, we can write 

d = 2 : p(x = x c ) = 1/F(0), i.e. p(x c , T) = const = p c . (36) 
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T 



FIG. 5: (ni'S) Temperature dependence of the resistivity p(T) at various values of x with the horizontal separatrix x = x c . (b) 
Reduction of the all the p(x, T) curves to two curves l)34[| of the dependence of the resistivity p scaling variable Il3"ill 

This implication (a horizontal separatrix in the set of temperature dependences at various values of x) is schematically 
shown in Fig. UK- We note that it is valid only under assumption (|25|) . which indicates that the correlation length £ 
is temperature independent. In general, p(x c , T) may have a finite slope at the point T = (cf. below, section 5.3 
and Fig. 13). 

Usually, with Sx replaced by the modulus \5x\, the p(u) curve is represented in the form of two branches as a 
function of In u (Fig. 5b) . The values of zv are then chosen such that these curves fit all the experimental points 
obtained at various temperatures. 

3. FLOW DIAGRAMS FOR METAL INSULATOR TRANSITIONS 

So far, we have not specified the type of a quantum phase transition. Hereafter, we speak about transitions related 
to a change in the electron localization, which, in turn, is caused by the degree of disorder in the system. The general 
theory of quantum phase transitions initially assumes that a control parameter affects the interparticle interaction and 
that a disorder is not more then a perturbative factor in the initial scheme of quantum phase transitions. Therefore, 
the applicability of this scheme to metal-insulator transitions, in which the degree of disorder is the main factor and 
the interaction is a secondary factor, is not obvious a priori. 

The fundamental difference between an insulator and a metal is that the electronic states at the Fermi level are 
localized in an insulator and delocalized in a metal. If an insulator is transformed into a metal due to a change in a 
certain parameter, the properties of the wave functions at the Fermi level change. The main physical property that is 
radically different in materials of these two types is the conductivity, i.e., the possibility of carrying an electric current 
at an arbitrarily weak electric field. This gives a 'yes-no' type signature: the conductivity is either zero, a = 0, or 
nonzero, though it may be arbitrary small. However, at a finite temperature T ^ 0, an insulator also carries a current 
owing to hopping conductivity. Therefore, the definition of an insulator given above is only related to the temperature 
T = 0, and the concept of a metal-insulator transition makes sense only at T = 0. 

That the conductivity a is not a function of the state (because it is realized only under nonequilibrium conditions) is 
an additional argument against the applicability of the general theory. However, Thouless Q noted that the transport 
properties can be used to characterize the decrease in the wave function of an electron placed at the center of an L d 
cube in a d-dimensional space when it moves toward the edges of the cube. Therefore, in analyzing the behavior of 
conductance or conductivity near the transition, we actually follow the evolution of wave functions. 

Thus, we have every reason to apply the general theory of quantum phase transitions to metal-insulator transitions. 
However, the first successful version of a theoretical description of metal-insulator quantum transitions |Q] was based 
on the renormalization group theory borrowed from quantum field theory. The essence of its conclusions, related 
to systems of noninteracting electrons in a random potential, is illustrated in Fig. El which shows the logarithmic 
derivative of the conductance of a sample with respect to the size L, 

p_ d\ny _L dy 
d In L y dL 

at the temperature T=0 as a function of the conductance y, 



(37) 
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FIG. 6: Universal f3(lny) functions for different dimentionalities 0|. 



The PdO^y) curves describe the universal laws of a change in the conductance of the system when its sizes change. 
The behavior of the system of noninteracting electrons in any material is described by a part of the curve (3d for the 
correspon ding dimension d. The interpretation and method of using these curves are described in detail in review 
and book [id ]. 

For metal-insulator quantum phase transitions, the theory 0] plays the role of an existence theorem. The shape of 
the (3d (In y) curves and their position in the (j/,/3) plane indicate that these transitions are absent in one- and two- 
dimensional systems of noninteracting electrons and that in three-dimensional systems, such a transition can occur, 
with the conductivity of the material changing continuously during this transition. This means that sufficiently large 
one- and two-dimensional samples should be always insulators at absolute zero whereas a three-dimensional material 
can be either an insulator or a conductor. Below, we discuss this feature in detail. Here, we note that the /3rf(lny) 
curves are called differential flow diagrams or Gell-Mann-Low curves. 

On the one hand, the theory in [l], 2, 0] is more specialized than that in 0] in some respects, because it only 
describes the vicinity of a quantum transition. On the other hand, the theory in [J, H, 0| is more universal, because 
its applicability is not limited by either interaction or a strong magnetic field. Therefore, it is of interest to compare 
the conclusions of these theories for the systems to which they both can be applied. 



4. THREE-DIMENSIONAL ELECTRON GAS 



The possibility of a phase transition in a three-dimensional material follows from the fact that the /^(lny) curve 
intersects the abscissa axis (3 = 0. A certain point on the /^(lny) curve (an image point) corresponds to the transport 
properties of a sample of size L made from a certain material (Fig. [7H). If this point is located in the lower half-plane 
(3 < the material is an insulator, and the image point shifts to the left as L increases; as a result, the conductance of 
the sample becomes exponentially small. If this point is located in the upper half-plane (3 > the image point shifts 
to the right; as a result, the material is a metal, and its conductance y increases with the size L. 

The curves in Figs [6] and EH are scarcely adapted for a direct comparison with experimental data because of a 
special choice of their coordinate axes. This disadvantage can be partly corrected by integrating the equation 



^=dlnL. (38) 

for three-dimensional (3D) systems. Since the left-hand side of differential equation f38|) becomes infinite at the 
point y = y c , the curves corresponding to different integration constants decompose into two families: one of them 
corresponds to an insulator region and the other to a metallic region. To demonstrate this by the simplest way, we 
restrict ourselves to the immediate vicinity of point y c , in which the /^(lny) curve can be approximated by a straight 
line, 

— — = s In — 39 

a In L y c 
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FIG. 7: Noninteracting 3D electron gas: (a) differential scaling diagram from [7| (see also Fig. [6}; the intersection of the curve 
/33 (In y) with the abscissa axis /3 = means existence of a quantum phase transition; (b) flow diagram for a noninteracting 3D 
electron gas in the vicinity of the transition point; (c) two universal curves of the dependence of the conductance of the 3D 
system on its size L obtained by scaling 



where s is the slope of the line with respect to the abscissa axis. Correspondingly, Eqn l|38|) can be replaced by linear 
equation lj39|) . The general solution of linear differential equation l(39|) is given by 



]n*=(f)V*, (40) 

where A plays the role of the initial condition that fixes the initial point in the ^(lny) curve; for example, it can be 
the length k^ 1 determined by the electron concentration. On the metallic side, fcpT 1 is equal to the minimum free 
path length l m i n . The conductance y\ corresponds to the point A. 

Figure [7b shows particular solutions of Eqn l|39p . If y\ > y c , we have lny\/y c > 0, the solution is located in the 
right quadrant, and the conductance y increases with L (metal). If y\ < y c , we have \ny\/y c < 0, the solution is 
located in the left quadrant, and the conductance y decreases as L increases (insulator). The set of curves in Fig. [7p 
is called a flow diagram. The arrows in this diagram show the direction of the image point motion when the system 
size increases. The flow lines corresponding to the particular solutions of Eqn lj39|) fill both upper quadrants of the 
(In yx/yc, L s ) plane. The boundary between them, the \ny\/y c — axis, is called the separatrix; in Fig. [7b, it is 
indicated by a dashed line. 

We extend straight line (f3Q|) in Fig. [7K to its intersection with the asymptote (3 = 1; as a result, we approximate 
the 03(y) curve in the upper half-plane > by a broken line consisting of segments of two straight lines. The size 
L which brings the image point to the intersection point is called the correlation length £. From Eqn l[40|). we have 

2A X 



£ = A(sln^) . (41) 

Vc 



With Eqn (|4Tj) . we can rewrite general solution lj40"l) as 



lnf=ff] (42) 



y c \£ 

normalize the size L by the length £ (which is specific for every material) , and reduce each of the two families in Fig. 
to one scaling curve (Fig. Et). 
We make several remarks. 

First, Eqns lf39f - lf42"j) involve a parameter s, which is unknown a priori. However, in the case of a noninteracting 
3D-gas, this parameter is either equal to unity or very close to it |9|. Because s enters all formulas as a factor or 
exponent, this parameter can be dropped. 

Second, the flow lines in Fig. ]7b are straight only because we integrated the Gell-Mann-Low equation only in the 
small vicinity of the transition point. If we integrated the ^(lny) function over a wide range of its argument Iny 
using some model representation, we would obviously obtain curved flow lines. 

The third note concerns the role of temperature. The parameter that specifies system motion along the flow lines 
indicated by arrows in Fig. [7b is the size. If the temperature is taken to be zero (as was assumed until now), this 
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parameter is the sample size L. However, we can initially suppose that the system is large, as this is done in the 
theory of quantum phase transitions [l|, 2, 0]. I n this case, the limiting size is considered to be the dephasing length 
L v , i.e., the size at which quantum coherence in an electron system is realized. The image point then moves along 
flow lines in the direction indicated by arrows as the temperature decreases. Because L v depends on the temperature 
and becomes infinite as T — > 0, this procedure introduces the temperature into the set of physical quantities related 
to the flow diagram. However, we do not study this relation; based on the 'existence theorem' of a quantum transition 
in a 3D space, we turn to the (x, T) phase plane to construct the critical region (see Fig. [8] and compare it with Fig. 

n. 

Let small values of the control parameter x correspond to metallic states. Then, at small x on the abscissa axis, 
the Drude formula a = 03 = ne 2 l/hk-p holds, and, at small x and finite T, the quantum correction 



a = a 3 + —L-\ L ee = y/hD/T. (43) 

is added to 03 (here D is the diffusion coefficient). The diffusion length in Eqn lj43j) is taken to be the length L ee 
determined by the Aronov-Altshuler effect [lH, i.e., by the electron-electron interaction. In other words, we assume 
that dephasing is caused by internal processes that occur in the electron gas in the time 

Tee = h/T, (44) 

without external impacts, such as the electron-phonon interaction (see analogous discussion of the physical meaning 
of L v in review [2], p. 324). 

At any control parameter x, another special point x m (the Mott limit) exists on the right of the transition point 
x c . At this special point, we have 

e 2 e 2 1 

a(x m ,T = 0) = o-Mott = -r-kp = —-. (45) 

In the range between x c and x m , the Drude formula is invalid: describing conductivity with this formula implies 
introducing a free path length / that is shorter than the de Broglie wavelength. In this range at T — 0, a is expressed 
through the correlation length £; this corresponds to the general scheme of the theory of quantum phase transitions. 
As a result, the conductivity along the abscissa axis is expressed as 

X ^ Xc 5 

17 [ T = in = ( — ) J L/i x c ^x > x m , (46) 

To match the last two expressions at x = x m , it suffices to set £(x m ) = I. 

Thus, to describe conductivity in the metallic region, we introduced two parameters, £ and L cc , which have the 
dimension of length and the properties required for the critical region (they are mutually independent and diverge at 
the point x c ). Following [l3|, we construct the following interpolation function in the critical region: 
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FIG. 8: Vicinity of a metal-insulator transition in a noninteracting 3D electron gas in the (x, T) phase plane 
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FIG. 9: (a) Temperature dependence of the conductivity of Ge:As samples with various level of doping in the region of a 
metal-insulator transition (from [3); the critical concentration is determined by the extrapolation of the experimental data 
to T = 0. b) Temperature dependence of the conductivity of a GaAs sample in various magnetic fields in the region of a 
metal-insulator transition (from [IB]); the critical field determined by extrapolation is 9.78 T. 



Expression (|47| matches Eqn (@3|) on the straight line x — x m and gives the correct values of conductivity in the 
gment at T = 0. The second form of the function demonstrates that function lf47j) satisfies general 

relation lj34|) . 

We now move from right to left along the line T = const (Fig. 01 line aa). As long as the quantum correction is 
relatively small, electron diffusion occurs via scattering by impurities; i.e., it is controlled by the first term in Eqn (|43|) . 
Therefore, the diffusion coefficient D in L ee is temperature independent. However, when we enter the critical region, 
03 transforms into (e 2 /fi.)£ -1 and begins to decrease rapidly. Under these conditions, D ceases to be a constant: 
diffusion is likely to be caused by the electromagnetic-field fluctuations that determine L cc ; as a result, this diffusion 
becomes temperature independent. We can then write a set of equations for the cx(T) and D(T) functions with the 
Einstein relation 




(48) 

used as a second equation. Here g-p is the density of states at the Fermi level. 

We eliminate D from two equations (|48|) . use that l/£ <C l/L ee near the transition, and obtain the temperature 
dependence of the conductivity on the right-hand side of the critical region l3, 13], 



aiT) = T G + {T9f)1/3 ) + /3Tl/3 ' (49) 

This means that inside the critical region, L cc is given by 

L ee = {Tgp)- 1 / 3 . (50) 



rather than being determined by Eqn (|43|) . Exactly at the transition, we have £ = 00 and a = 0. 



Temperature dependence (|49|) was repeatedly detected in experiments performed on a variety of materials [14j, [15J, 
16 , [rzj. Figure [9] shows two examples in which the control parameter is represented by the electron concentration and 
magnetic field. 

Using condition l(32|) and Eqn l|49p . we can write the relation 

T oc G/ff 8 )- 1 . (51) 
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for the right boundary of the critical region. If the relation £ oc ((Sx) _1 holds (i.e., if v = 1, as is assumed in Fig. [8]), 
the boundary of the critical region is represented by a cubic parabola T oc (Sx) 3 . In the general case, we have v ^ 1 
and 

T oc (5xf»/g F . (52) 

A comparison with Eqn lj33j) demonstrates that in the metal-insulator transition in a 3D system of noninteracting 
electrons, the dynamic critical index is 

2 = 3. (53) 



Below curve |52j) in the x c > x > x m segment, the conductivity is also described by Eqn (|47|) : however, L cc enters 
Eqn f47j) in the initial form f43|) with the diffusion coefficient -D = const. Therefore, the temperature dependence of 
the conductivity in this region should have the form cr(T) = a + fiT 1 / 2 . 

Up to this point, we have discussed the right-hand side of the phase diagram. The left-hand side of the phase 
diagram, i.e., the insulator region (where hopping conductivity takes place), also has two characteristic lengths. First, 
there is the decay length £ of the localized-state wave functions, ip oc exp (— r/£). Far from the transition, £ decreases 
to the Bohr radius = nh 2 /(m*e 2 ) (where k is the dielectric constant and m* is the effective electron mass); at 
the transition, it diverges because the electrons become delocalized. Second, there is the average hopping distance r. 
If the hopping conductivity is described by the Mott law, the average hopping distance is [10| 

r={£,/g B T) 1 ' i . (54) 

The lengths r and £ cannot be used as two independent lengths in the critical region because they are connected by 
relation lj54j) . However, using L ee determined by Eqn |50| . we can rewrite Eqn (|54| as 

r = m e ) lf \ (55) 

We take into account that expression (|50|l for L cc does not contain the kinetic characteristics of the electron gas and 
consider L cc the dephasing length over the entire critical region, including its left-hand side, above the x; > x > x c 
segment on the abscissa axis. Equation (|52"1) then determines both the right and left boundaries of the critical region. 
Equations 1|54[1 and (f55|) demonstrate that r = £ = L cc at the left boundary. 

The conclusion that L cc can be considered to be the dephasing length over the entire critical region is supported 
experimentally: as can be seen from Fig. [9l the a(T) temperature dependence straighten in terms of the (T 1 / 3 ,^) 
axes, not only on the right-hand side of the critical region but also at values of the control parameter x < x c (e.g., 
at the electron concentration n = 3 x 10 17 cm~ 3 in Fig. [9]). However, the free term a in Eqn f49|) becomes negative. 
This means that we should either suppose that correlation length £ is negative in the insulator region or (which is 
formally preferred but is essentially the same) should replace interpolation formula |49|) on the left-hand side of the 
critical region by the formula 

As follows from Eqns l|5"5"Tl and <[56|l . we have a — along the left boundary of the critical region, which means that 
the conductivity along this boundary is determined up to an exponentially small hopping conductivity. 

The hopping conduction mechanism is still operative below the lower boundary of the critical region over the 
x\ > x > x c segment, and the wavefunction decay length is £ 3> ae rather than ae [13]. Therefore, ae does not enter 
the expression for the hopping conductivity, and a is expressed in terms of the correlation length £ (see Fig. [8]): 

a oc exp(— r/£). 

The (x,T) phase diagram in Fig. [8] accumulates the results of the long-term experimental studies of the low- 
temperature transport properties of conducting systems, namely, the quantum corrections to metallic conductivity, 
the evolution of electronic spectra during metal-insulator transitions, the temperature dependence of conductivity 
in the vicinity of the transitions, and hopping conductivity. In essence, this phase diagram was plotted irrespective 
of the theory of quantum phase transitions |lj, 2, 3] in order to reveal the compatibility of all experimental data. 
Nevertheless, the considerations given above demonstrate that this diagram is absolutely adequate for the concepts 
following from this theory. 
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5. TWO-DIMENSIONAL ELECTRON GAS. 
5.1 Gas of noninteracting electrons 

We now pass to 2D systems. According to the theory by Abrahams et al. [7], a 2D system (d = 2) of noninteracting 
electrons is always an insulator in the sense that localization should inevitably occur in a sufficiently large sample, 
L > £, at a sufficiently low temperature, T < T^, and an arbitrarily small disorder. This statement follows from the 
fact that flow line d = 2 in Fig. [6] asymptotically approaches the (3 = axis and does not intersect it (this line is shown 
in Fig. [TOh). However, in weakly disordered films, the correlation length £, which bounds the size L below, can be 
unrealistically large, and the temperature of the crossover from the region with a logarithmic quantum correction 
to the conductivity to the region that is characterized by an exponential temperature dependence is, in contrast, too 
low. These films are called metallic films. 

Estimates of £ and XJr can be obtained from the assumption jitj that the logarithmically diverging quantum 
correction Ac in the conductivity 

a = cr 2 + Act = (e 2 /h){k F l) - {e 2 /K) \n(L ee /l) (57) 

is of the order of the classical conductivity cr 2 , and hence a ~ and 

HL ee /l) « k F l. (58) 

The diffusion length L cc over which electron dephasing occurs is determined by Eqn (j43| . The value of L ee 
determined from Eqn (|55)l is defined as the correlation length £, 

£ = Zexp(M), (59) 

and the temperature 

T 6 = Dh/e, (60) 

determined from the relation L cc = £ using Eqn lj43j) is called the crossover temperature. Of course, the conductivity 
does not vanish at T — Tg; however, the theory of weak localization is obviously not valid at this temperature, and 
the resistivity should begin to increase exponentially with decreasing the temperature at T < in samples of size 

We now plot function ([60)) and, for convenience, direct the l/£ axis to the left (see Fig. [10]). As a result, this 
diagram can be conveniently compared with the diagram for a 3D system shown in Fig. El We add an axis and lay 
off the control parameter x as abscissa; this parameter characterizes the degree of disorder, with x = corresponding 
to an ideal system in which a disorder is absent. 

It is easily seen that the T(l/£) curve in Fig. [10] represents the left-hand side of the phase diagram in the vicinity 
of a metal-insulator transition in a 3D material (cf. Fig. [8|) whose phase transition point is located at the edge of 




FIG. 10: Noninteracting 2D electron gas: (a) differential flow lines for systems with dimension d = 2 taken from the scaling 
diagram in [3] (see also Fig. EJ) and d — 2 + e (see text), and (b) crossover from the logarithmic temperature dependence of 
conductivity to its exponential dependence, which can be treated as the boundary of the critical region of a virtual quantum 
transition (see text) 
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the diagram, at the origin (T = x = 0). Curve f60|) is then the left boundary of the critical region, and the dynamic 
critical index is 

z = 2. (61) 

Using Eqns lf53|) and <f6T]) . one can conclude that for metal-insulator transitions in systems of noninteracting 
electrons, the dynamic critical index is equal to the dimension, z — d. 

Thus, systems with dimension d = 2 turn out to be boundary systems: a quantum transition is still present in 
the {x, T ) phase plane but is shifted toward its corner, to the unreachable point x = 0. The boundary properties of 
2D systems can also be found from the flow diagram in Fig. [6J We imagine that the dimension d is a continuous 
parameter and can take not only integer values. Straight lines (3 — d — 2 represent the asymptotics of the flow lines at 
high values of y. Therefore, the flow lines of a system with dimension d = 2 + e > 2 inevitably cross the abscissa axis 
(3 = 0, and such systems have a metal-insulator transition (Fig. flOa). As e — > 0, the transition point shifts toward 
high conductance and goes to infinity. 

This interpretation of the curve in Fig. [TO] implies that the domain over the T cx (l/£) 2 parabola is the critical 
region of the quantum transition. For 2D systems, only one scaling variable u is retained in Eqn f34|) written for the 
conductivity in the critical region; it is equal to the ratio of two characteristic lengths, 

o = F(u) = F{L V I£). (62) 

In this region, however, the conductivity is typically expressed as the difference between the classical conductivity 
and the quantum correction, 

a = a 2 - Act = ^ (k F l - In ^\ . (63) 

To resolve this apparent discrepancy, let us substitute the expression I = £ exp (— kpl) from Eqn 1(59]) in the argument 
of the logarithm in Eqn lj63j) . move the exponent from the logarithm, and obtain 

The classical conductivity o-i cancels and the remaining part depends only on the scaling variable, as it should be in 
the critical region. Thus, in this regard, our (x,T) diagram also satisfies the requirements of the theory of quantum 
phase transitions. 

Expression f64|) holds not in the entire 'quasi-critical' region T > Tj; moreover, it diverges on the x = axis. 
However, near this axis, the elastic mean free path I — > oo and the standard expression for L cc lose their meaning 
because the necessary condition t cc 3> t is violated. The boundary of the part of the region where Eqn (|64| is invalid, 
t oc = r, is shown in Fig. H0]by a dashed straight line plotted under the assumption that r cx a; -1 . 

As is shown in Section 5.3, the introduction of interaction can lead to the appearance of a metal-insulator transition 
in a 2D system. In terms of the phase diagram shown in Fig. [TOj this means that the phase transition point shifts 
from the origin to a point x c / on the abscissa axis due to a certain cause, and a boundary additionally appears 
between the critical and metallic regions. Standard expression ([63| . which was written for the conductivity in the 
metallic region and was transformed into Eqn ([64]) . was already used for the description of the conductivity in the 
critical region. This means that the expression for the conductivity in the metallic region should be radically different 
and that we should expect the appearance of a marginal metal instead of a Fermi liquid. The nature of the electronic 
states in this hyp othetic metal should be rather peculiar, since electrons are assumed to be localized when interaction 
is turned off jilf. 

5.2. Spin— orbit interaction 

The boundary position of 2D systems makes them sensitive to various types of interaction, e.g., spin-orbit interaction 
or electron-electron interaction, which can cause a phase transition. 

We first consider the spin-orbit interaction. This case is convenient because a flow diagram can be plotted using the 
(y, /3) axes of the initial diagram shown in Fig. [6l At large y, the initial flow line /?2 (y) deviates down from its right 
asymptote j3 = because of weak localization, which results in a decrease in the conductivity (Fig. fTOk). But the 
spin-orbit interaction changes the sign of the quantum correction, i.e., changes weak localization into antilocalization. 
Therefore, the sign of the derivative on the right-hand side of the /^(lny) curve should change: the curve deviates 
upward from the asymptote (3 = 0, goes to the upper half-plane (3 > 0, and, hence, inevitably crosses the abscissa 
axis (3 = and reaches the left asymptote. 
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Antilocalization was studied in detail both experimentally ji^l and theoretically j23|. In the 2D case (d = 2), the 
quantum correction to the conductivity is given by 

A«»-<™/f(i^-i). («) 

T 

where r is the time between elastic collisions, r so > r is the time between spin-flip collisions, and r oc = H/T is the 
dephasing time. The motion of the image point in the flow diagram now depends on two parameters, the conductance 
y and time r so . For various values of r so , Fig. fTTI shows the family of flow lines located between two envelope curves. 
The lower curve is fa (In y) from the diagrams in Figs[6]and[T0l It can be obtained from integral f65|) if we set r so = oo 
in the integrand, and hence the parenthesis in the integrand becomes equal to unity. The upper curve has the same 
asymptotics; however, its construction implies that r so = r. Then, we have t » r so over the major portion of the 
integration range, and the parenthesis is considered to be —1/2. 



p 







— *\ — 



FIG. 11: Flow diagram for a 2D electron gas with spin-orbit interaction. 
We now take a sample of size L of a material with an intermediate value of t so , 

T C T so C OO. 

The size L bounds the diffusion time t of an electron in a 2D sample until the collision with boundary by the quantity 
TL,t < tl ~ t(L/1) 2 . Therefore, to describe weak localization in a sample of size L, we must replace the upper 
integration limit in Eqn lj65ll with tl, 

Tee -> T L = t(L/1) 2 . 

Let L be first very small (L 

gtrsiml) and the temperature T = 0. Diffusion processes then have no time to develop; electron interference is 
virtually absent; the conductivity is equal to its classical value; and the image point is on the right on the axis = 0. 
As the size L increases to 

L so = W t so/t (66) 

spin-flip collisions are insignificant and correction lj65j) to the conductivity is negative. The image point moves to 
the left along the lower envelope curve, as in Fig. [TO] in the absence of the spin-orbit interaction. When the size L 
becomes larger than L so , integral J65|) and A02 change their sign and the image point moves to the upper envelope 
curve. If this passage occurs at y < y c , the film becomes an insulator at L — > 00 in any case. But if the spin-orbit 
interaction is strong, r so is small and the image point reaches the upper envelope curve at y > y c . Then, the image 
point continues to move along the upper envelope curve toward large y, and the film retains its metallic properties as 
L — > 00. Thus, a metal-insulator transition could be observed in an experiment with the spin-orbit interaction used as 
a control parameter. So far, only the transformation of weak localization into antilocalization has been demonstrated 
this way [22|. 

5.3. Gas of interacting electrons 

We have to refine what is meant by the absence or presence of the electron-electron interaction. The interaction 
manifests itself differently in the properties of electron systems; for example, it determines the probability of electron- 
electron scattering. In the classical theory of metals, electron-electron scattering is considered not to contribute to 
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conductivity, since the total momentum of the electron system and the drift velocity remain the same. When quantum 
effects are taken into account, this statement becomes invalid because the conductance depends on the dephasing 
length and electron-electron scattering changes this length. Nevertheless, if the electron-electron interaction affects 
the conductance only through scattering, the electron system is considered noninteracting, because the interaction 
affects the conductance as an external action, e.g., scattering by phonons. 

Electron-electron scattering is not the only channel of the effect of the interaction on the conductance. The 
interaction is thought to be not the bare Coulomb interaction but the interaction between 'dressed' quasiparticles, 
i.e., the screened interaction that depends on the electron density, the diffusion coefficient of electrons, and (under 
certain conditions) the sample size or the quantum coherence length [24 1. 

This interaction determines the structure of quantum levels and the ground-state energy, affecting the competition 
between phases in the vicinity of a phase transition. Electrons are mainly scattered by a screened impurity potential, 
whose properties depend on the effective electron-electron interaction. This interaction, in turn, depends on the 
diffusion coefficient and dephasing length, which results in the complex dependence of the effective interaction on all 
external parameters. 

The case of spin-orbit interaction considered in Section 5.2 represents an example of interaction renormalization 
as the sample size L changes: the spin-orbit interaction is actually turned on only when L becomes larger than the 
value determined from Eqn (|66|) . 

As another example, we analyze the picture of possible states in the model with a multivalley electron spectrum 
developed by Punnoose and Finkel'stein In this model, a change in L causes changes in both the conductance y 
and interaction 9. Therefore, an equation for 9 is added to Eqn l(37|). The flow diagram is a result of the solution 
of these two equations. Figure fl2l shows part of this diagram calculated by the authors of in slightly different 
coordinate axes. The variables were changed to facilitate a comparison of this diagram with those shown in Figsni'S6 
and 10, although the quantitative information contained in the initial diagram in (25J is lost. The abscissa of the flow 
diagram shown in Fig. [12] is the conductance. Thus, the abscissa of all the diagrams in Figs [6] and [10] is the same. 
The ordinate 9 of the flow diagram in Fig. [12] reflects the effective interaction. A noninteracting electron gas, which 
was discussed in Section 5.1, corresponds to the straight line 9 = 0. As in all previous diagrams (Figs I MTOll . the size 
is the parameter that determines system motion along the flow lines indicated by arrows. This size is given by the 
sample size L if T — or by the dephasing length L v (or L cc ), i.e., the maximum size in which quantum coherence 
is retained in an electronic system. 

The configurations of the flow lines in Fig. [12] clearly display the specific feature of the interaction 9 that was 
discussed at the beginning of this section: as the size L or L v changes, the interaction in the system changes (is 
renormalized), and this change is different in different flow lines. The diagram in Fig. [T2] is a two-parameter diagram: 
two parameters are required to specify the flow line of the image point. The flow lines (trajectories) in this diagram 
occupy the entire half-plane. 

We assume that the position of a point in a flow line is determined by L v , i.e., by the temperature. We assume that 
we are at point A in a flow line, that is metallic because it recedes to the region of high conductivity at T — > 0. We 
fix the interaction 9 and vary the conductance y; that is, we vary the degree of disorder. This process is shown by a 
dashed straight line in the diagram. When moving along this line, we can cross the separatrix and reach point B in a 
flow line that describes an insulator, because it tends to the point y = as T — > 0. To specify the position of a point 
in line AB, we can use a single parameter, which is called a control parameter. At another control parameter, the 




FIG. 12: Part of a flow diagram plotted for a 2D electron gas with interaction [25|] 
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angle of intersection of the separatrix can be different. For example, we can initiate the crossover from a metallic to 
a nonmetallic flow trajectory by changing the electron concentration and, thus, simultaneously changing the effective 
interaction and the conductance. 

By changing the control parameter, we can pass from one flow trajectory to another, and, by changing the temper- 
ature (size), we can move along a flow trajectory. But because interaction 8 can change when any of two parameters 
changes, factorization is absent; that is, we cannot suppose that the length £ depends only on Sx and the length L v 
depends only on temperature. The consequence of this 'mixing' of variables is clearly visible in the flow diagram in 
Fig. Q21 As we move along the separatrix toward the quantum critical point (QCP), the conductance decreases and 
tends to y c - Hence, we can draw the following important conclusion, which is qualitatively shown in Fig. [131 the 
separatrix in the set of temperature dependence is not a horizontal line, as in Fig. [H 



y, 




FIG. 13: Qualitative scheme for the evolution of the conductance y(T) curves with a control parameter in the model proposed 
in 1251 



An indication of the presence of a metal-insulator transition in a 2D gas was first obtained in the inversion layer of 
a field-effect transistor on a Si surface |2f|. The presence of the transition was questioned for a long time, because is 
was in conflict with the concepts formulated in [7|| and because the transition was not reproduced in other materials. 
However, the uniqueness of silicon was found to be related to a high electron mobility, which allows performing 
experiments at a very low electron density, where the electron-electron interaction is especially important. With 
this fact, it was possible to interpret the experimental data using not the theory in Q, which was developed for 
noninteracting electrons, but in accordance with the model in [25j. As a result, the presence of a metal-insulator 
transition was s upp orted and a finite slope of the separatrix, which was predicted in [25|], was obtained (see Fig. [T4r 
borrowed from [27j. 

As is seen from the flow diagram in Fig. 12, the finite slope of the separatrix in the set of the temperature dependence 
of a or R of a system of 2D electrons is controlled by the angle at which the separatrix in the flow diagram approaches 
the QCP. If the tangent to the separatrix is normal to the abscissa axis at the QCP due to any specific reason, the 
separatrix in the set of temperature dependence has a zero derivative at T = 0. Thus, the horizontal position of the 
separatrix in the set of the temperature dependence of conductivity in the vicinity of the quantum phase transition 
results from the symmetry of the flow diagram of a certain system and is not an inherent property of all 2D systems. 
Conversely, a finite slope of the separatrix is not an indispensable consequence of a two-parameter flow diagram. 

An inclined separatrix complicates scaling, i.e., the reduction of measurements performed along different flow lines 
to one universal curve by changing the scales. Nevertheless, the scaling of the resistance R{T) data is possible. Figure 
[T4b shows the scaling carried out in the metallic region of the transition displayed in Fig. [T4h . The three lower curves 
in Fig. fl4k are replotted in the coordinates p/p max (instead of p ) and /9 ma xhi (T/T max ) (instead of T; here, p max is 
given in dimensionless units) , and the values p max and T max correspond to the position of the maximum in each of 
the experimental curves, which are seen to merge into one curve and to coincide with the theoretical curve. This last 
curve was plotted using the calculations in [H] and the interaction parameter 9 determined in the same sample from 
the magnetoresistance data obtained in a field parallel to the 2D plane. 

An analogous problem of an inclined separatrix is also often encountered during the experimental processing of 
the R x {T ) temperature dependence (x is a control parameter) in the vicinity of superconductor-insulator quantum 
phase transitions. In Ref. [28], a procedure was proposed for the 'correction' of the curves via the introduction of the 
additional linear term 



R X {T) -» R X {T) - aT, a 



dT 



(67) 



T=0 



into each of them to make the separatrix horizontal and for performing subsequent standard scaling for 2D systems 
using scaling variable ([35| . The idea of this procedure is to compensate for the slope of the separatrix in a flow 
diagram. However, the correctness of this procedure has not been proved theoretically. 
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FIG. 14: a) Temperature dependence of the resistivity of a 2D gas in a field-effect silicon transistor in a concentration range 
containing the metal-insulator transition [27l |. b) Scaling of the three lower curves, which correspond to the concentrations 
9.87 (□), 9.58 (O), and 9.14-10 10 cm" 2 (+). Theoretical results are shown by a solid line, [13] 




FIG. 15: Critical region in the flow diagram of a 2D electron gas with interaction (indicated by a dotted curve) 



Formally, the (x,T) phase plane also has a meaning for a two-parameter flow diagram. However, the critical 
vicinity of the transition can also be directly plotted in the diagram (see Fig. 1X5]) . It should be noted that the 
two-parameter scaling becomes one-parameter in the region adjacent to the separatrix inside the critical region, where 
flow trajectories are parallel. Indeed, as noted above, one parameter is sufficient to specify the position of a point on 
line AB, and this parameter does not change when line AB shifts parallel to the separatrix, e.g., to position AiBi, 
since the flow lines are parallel to each other. As we move along the separatrix toward the quantum critical point, 
near-separatrix trajectories turn aside alternately, and the strip in which the flow trajectories are parallel, as well as 
the critical region, narrows. 



6. QUANTUM TRANSITIONS BETWEEN THE DIFFERENT STATES OF A HALL LIQUID 

States in the plateaus of the quantum Hall effect are the specific phase states of a 2D electron gas with special 
transport properties described by the longitudinal a xx and transverse a xy conductivities 

a xx -^0 at T™>0, a xy = i{e 2 /2Trh), i = 0, 1, 2, 3... . (68) 

Such phase states of a 2D electron gas are quantum Hall liquids with different quantum Hall numbers i, which are 
determined by the values of the Hall conductivity a xy , Eqn f68j) . in the plateaus, 

i = a xy /(e 2 /2irh), i= 1,2,3.... (69) 

The transitions from one plateau to another induced by a change in the magnetic field or the electron concentration 
are clearly visible in experimental curves. As can be seen from Fig. [TBI a xy jumps are accompanied by narrow a xx 
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peaks. These jumps are quantum phase transitions, and they should fit both theoretical versions compared in Section 
5. 
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FIG. 16: The magnetoresistance R xx and the Hall resistivity p xy versus the magnetic field B in a GaAs-ALGai-x As het- 
erostructure at T = 8 mK [29]. The electron density in the 2D layer is 3.7 ■ 10 11 cm -2 , and the mobility is fi = 4.1 • 10 4 cm 2 /V-s. 



Figure [TBI only displays the plateaus of the integer quantum Hall effect, which is considered below. The integer 
quantum Hall effect can also be realized in a noninteracting electron gas; therefore, it can be described without regard 
for interaction. 

To construct the flow diagram of a 2D system of noninteracting electrons in a strong magnetic field, we need two 
conductance components that are equivalent to conductivity components a xx 
transforms into the set of two equations 



and <j X y Correspondingly, Eqn (j37l 



d In <t xx 
d\nL 

d In <7 xy 
d\nL 



~ fl faxx i @xy ) i 



— f2(&xx,0- xy ), 



(70) 



When eliminating the variable L from two equations (|70|) . we find a relation between <r xx and <r xy , which can be 
displayed as curves in the (a xy , a xx ) plane [30]. These curves make up a flow diagram for a 2D system of noninteracting 
electrons in a strong magnetic field (Fig. [T7j) . The flow lines in this diagram are separated by separatrices periodically 
repeated along the a xy axis. This diagram is again a two-parameter diagram, but due to a strong magnetic field 
rather than to interaction. As in the previous flow diagrams, we can move along arrows in flow trajectories by either 
increasing the sample size L or decreasing the temperature at large L, i.e., increasing L v . 

In the region below points Q, the separatrix is split such that two equivalent points located at the same height 
(Bi and B2) appear in it. The flow lines inside the semicircles in Fig. 17 are omitted because they are thought to 
exist separately from the flow lines outside these semicircles: as the control parameter in the region below points Ci 
changes, the motion leads to a hop between points (Bi and B 2 ) and to a sharp increase in a xy . Strictly speaking, 
the transfer should occur along a line located outside the semicircle and bypassing the point C; above. However, for 
simplicity, it is indicated by a horizontal dashed line. 

The control parameter in the quantum Hall effect regime is usually given by the electron concentration or the 
magnetic field. Their effect on the state of a real system depends on the random field of impurities and other defects, 
which transforms discrete Landau levels into minibands and specifies the energy structure and the character of wave 

1 /2 

functions in them. Because a magnetic field determines the magnetic length tb = (hc/\e\B) as a characteristic scale, 
we can speak about two limiting types of random potential, a potential with large-scale fluctuations of characteristic 
sizes C ^ r B and a short-range potential with (<Tfl. In the long-period potential model, an energy value e c exists 
near the center of each Landau miniband such that a delocalized electron wave function corresponds to this value. If 
the wave function is strictly delocalized only at £f = e c and the random field lifts the degeneracy of energy levels, a 
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FIG. 17: Flow diagram for a 2D electron gas in a strong magnetic field [30]. The coordinates are represented by the conductivity 
tensor components o xy and a xx in the dimensionless units e 2 /2-kH. Separatrices are indicated by dashed lines. Aj are stationary 
singular points, C% are unstable singular points, which are quantum transition points similar to QCP in Fig. 1121 The dashed line 
indicates the critical region near C2. Horizontal dashed line B1-B2 indicates the image point hopping as the control parameter 
changes. The right-hand side of the figure shows the hypothetical flow diagram that corresponds to a split phase transition 
and the appearance of a metallic phase 



smooth change in the electron concentration produces a jumplike transformation of the system from one phase state 
into another through an isolated energy state with a delocalized wave function at the Fermi level. This behavior is 
implied in the flow diagram in Fig. 17, where each isolated metallic state corresponds to a peculiar separatrix. 

The actual width Se of the energy range with delocalized wavefunctions depends on finer processes, such as tunneling 
between two semiclassical trajectories that are close to each other in the vicinity of the saddle point (magnetic 
breakdown). In essence, Se is the energy uncertainty of any delocalized state. Another source of increasing the Se 
range is the finiteness of the lengths L and L v . If the Se range is finite, a separatrix is split into two parallel lines 
and the phase transition is split into two transitions: a metallic state with a partly filled layer of extended states 
at the Fermi level should appear between two Hall-liquid states whose indices i differ by unity. The corresponding 
hypothetical flow diagram is shown on the right-hand side of Fig. 17. 

At first glance, it seems that experiment can distinguish between these two hypothetical possibilities. For definite- 
ness, we assume that the electron concentration n changes in experiment (considerations for a change in the magnetic 
field are similar). As the concentration changes, the Fermi level moves along the energy scale. When states at the 
Fermi level are delocalized, the a xx conductivity is finite, and the a xy conductivity is in the intermediate region 
between two plateaus. Therefore, the temperature dependence of the concentration range of the intermediate region 
(i.e., the a xx peak width, the da xy /dn derivative at the center of the intermediate region, etc.) extrapolated to T = 
must determine the energy range Se of the delocalized states. 

However, the experimental results were found to be ambiguous. On the one hand, the comprehensive experiments in 
[3l|,l32| give a finite value of Se. For example, Figfl8l shows the transition width measured from the p xx resistivity peak 
width of a 2D gas in the GaAs/AlGaAs heterojunction. The measured function is seen to be reliably extrapolated 
to Ai?o « 0.35 as T — > 0. However, as we see below, many experiments give the opposite result: as the temperature 
decreases, the transition width tends to zero (see, e.g., Fig. 120)1 . The factor that determines Se has not yet been 
revealed. In any case, there is no unique relation between Se and mobility. The statistical characteristics of the 
random potential, which are poorly controlled, are likely to play a key role. 

We now move to experiments that do not exhibit a finite energy layer with delocalized states. Although the flow 
diagram of i — » i + 1 transitions between different states of a Hall liquid has two parameters (see Fig. [17)), it is 
symmetric with respect to both the a xy — (i + l/2)(e 2 /2nh) and a xy — i(e 2 /2wh) axes. Therefore, we may use 
the version of the general theory of quantum phase transitions that is based on Eqns (|25)l . (|26)l and assumes that 
scaling formulas of type (|34)l are valid. As far as we consider only 2D electronic systems, all resistances have the same 
dimension [SI] and must have the form 

R U v — F uv (L(p/£) — F uv (^rpi/ zl , 1 > (71) 

in the vicinity of the transition. Here, the u and v subscripts stand for the x and y coordinates, F uv is an unknown 
function, and the argument of the arbitrary F function is written using the last form in Eqn (|3o)) . 
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FIG. 18: Temperature dependence of the peak width of the longitudinal resistivity p xx of a 2D electron gas in the GaAs/AlGaAs 
heterojunction during the 2^1 transition in the magnetic field about 4ni'S |32l | . The carrier mobility and concentration at 
T = 1.5 K are /i = 34000 cm 2 /Vs and n = 1.4 • 10 11 cm -2 . The inset shows p xx {B) curve recorded at 150 mK 



In contrast to the case of metal-insulator transitions, we analyze experimental results instead of calculating or 
predicting the values of v and z. As an example, Fig. [l9]shows the magnetic-field dependence at various temperatures 
of the longitudinal R xx and transverse R xy resistances of a Hall bar in a GaAs-based heterostructure [33j . 




5(T) 5(T) 

FIG. 19: Transverse R xy and longitudinal R xx resistances of an AUGai-a; As— Alo.33Gao.67As heterostructure, x = 0.85%, at 
various temperatures. The critical magnetic field of the 4-3 transition determined from the point of intersection of the R xy {T) 
curves is B c = 1.40 T, ^ 

According to Eqn lj7i~j) , the scaling variable 

u = 8x/T 1 ' zv (72) 

is identically zero at all temperatures in the case where the control parameter takes a critical value; correspondingly, 
we have 

R uv {x c ,T) = const. (73) 

Therefore, separatrix (73) must be horizontal, and all isotherms R uv (x,T — const) must intersect at one point x — x c . 
This is the first test of the applicability of Eqn (|7T]l . 

We first focus on the R xy (T) curves. As is seen from Fig. HH the R xy {T) isotherms obtained for an Al^Gai-^As- 
Alo.33Gao.67As heterostructure with x — 0.85% [33| do intersect at one point, B c = 1.40 T. Near the intersection 
point, all the R xy {T) curves can be expanded into a series and replaced by straight lines (dR xy /dB)B c (B — B c ). As 
the slopes of these lines are changed from (dR xy / '3B)b c to {dR xy /dB)B c /T K , where k = 1/zv, all the straight lines 
must merge into one line. The choice of the value of k at which the relation 



(dR xy /dB) Bc (T)/T K = const, 



(74) 
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FIG. 20: Determination of the critical index k = 1/zu for transition 4-3 using the R XX {T) peak width (a) and matching the 
slopes of the intersecting R X x{T) curves (b). The data were obtained from the sample used for Fig. 19 in a dilution refrigerator 
(triangles) and a liquid 3 He cryostat (circles) [33] 



holds is the second step in the application of the scaling procedure, and the possibility of this choice is the second 
condition for the applicability of the theory. To choose n, we plot (dR xy /dB) versus T on a log— log scale (Fig. l20j) . 

Formally, Eqn l(7T|) can be applied to both longitudinal (R xx ) an d transverse (R xy ) resistances. However, the peak 
height depends on the temperature; that is, the point of the maximum does not satisfy condition <[73j) . Leaving aside 
the reasons for this fact, we can use the longitudinal resistance data for scaling analysis by accounting for the integral 
property of the R XX (B) functions in the vicinity of the transition, namely, the peak half- width AB determined by a 
certain algorithm. In Fig. [2Ql AB is determined as the distance between the two maxima of the (dR xx /dB) derivative. 

As can be seen from Fig. [2Ql an analysis of both families of the functions gives the same value of the critical index 
k = 0.42, which is an additional argument for our scaling procedure. 

This value of the critical index k — 0.42 was repeatedly obtained in heterostructures made from various materials. 
However, contrary to expectations, this value is not universal: other values in the range from 0.2 to 0.8 were also 
detected in a number of experiments on various heterostructures. This scatter calls for explanations, since scaling 
relations and critical indices are usually universal. 

At T = 0, transition occurs when the condition ef — £ c is satisfied and therefore the difference Se = |ef — £c\ is 
the only 'internal' control parameter of the system, with the correlation length £ depending on it in accordance with 
a power law: 

foe (&)"*. (75) 

From this standpoint, the magnetic field B or the 2D electron concentration n, which depends on the gate voltage 
V g , are 'external' control parameters x, 

5x = \B C - B\ oc (5e) U2 or Sx = \n c - n\ oc \V gc - V g | oc {Se)" 2 . (76) 

In both cases, the exponent is the same. This is supported by the fact that the R UV (B) and R U v(V g ) experimental 
curves recorded using the same sample under equivalent conditions differ only in the scales on the abscissa axis. 
Eventually, we can write Eqn |25|) as 

€ cx {8xY {v = v x v 2 ), (77) 

where the relation between the control parameter 5x and the correlation length f consists of the following two links: 
f depends on the position of the Fermi level £f with respect to the delocalized level e c , and the difference ef — e c , 
in turn, depends on Sx. Correspondingly, according to Eqn |77| . the index v turns out to be the product of v\ and 
V2- The relation between f and 6s and the related index v\ are likely to be universal and the same for all transitions 
between different quantum Hall liquids. But the index vi is determined by the density of states g(e) in the vicinity 
of the energy e c and, hence, depends on the specific features of the random potential; this potential can be lon g- o r 
short-range, statistically symmetric or asymmetric with respect to the mean value, and so on. The authors of [33j], 
whose curves are used in this section, just studied the effect of specific features of the random potential on k. 

The AB(T) curve in Fig. l20l differs radically from the curve in FigHHJ in the latter case, a free term ABq was 
introduced for the experimental data to be approximated by a power function. Nevertheless, we can also perform 
scaling analysis of the experimental data in this case using the second hypothetical version of the flow diagram (see 
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Fig. fTTj) , For this interpretation, the presence of the ABq term means that in the 2 ABq range of the control parameter 
B, the image point in the flow diagram moves across the metallic-phase corridor and falls on a separatrix not at Bq 
corresponding to a maximum of a xx (B) and the derivative da xy /dB but at B + AB . This problem is discussed in 
detail in review Here, we only note that Fig. [18] actually contains this scaling analysis. By extrapolating the 
AB(T) dependence to T = 0, we can determine the delocalized-state layer width in units of magnetic field 2ABq and 
find that AB(T) ~ ABq depends linearly on T. This means that K = 1 in this experiment. The same value of k was 
obtained earlier in |31| . 



7. CONCLUSION 



All the considered cases of metal-insulator transitions were found to be adequately described by flow diagrams. 
The list of theoretical works that have successfully used this technique begins with work Q, where a theoretical 
model for noninteracting electrons in a zero magnetic field was developed. The last achievement in this field is the 
construction of a flow diagram for a 2D model system of interacting electrons and the demonstration of the possibility 
of a metal-insulator transition in this system . 

As regards the general theory of quantum phase transitions, a phase transition in a 3D system of noninteracting 
electrons demonstrates that, in principle, this theory can be used to describe localized-delocalized electron transitions. 
Neither conductance, which is used as a physical quantity specifying the state of the system, nor disorder, which is the 
main control parameter, are substantial obstacles for this theory. However, as usual, various particular cases require 
theoretical versions of various de gree s of complexity. For example, the version described in this review cannot be 
applied to the model proposed in [251 ] . 

The relative role and possibilities of both theoretical approaches are demonstrated when the integer quantum Hall 
effect is described. The flow diagram in Fig. [17] is very convenient for the discussion of the types of transitions that 
are possible in a system and for the formulation of questions to be experimentally checked. Many of these questions 
are still open. For example, it is unclear which of the versions of the flow diagram in Fig. [17] is realized in reality and 
whether the transition between the states with quantum indices i and i ± 1, 

% T± i ± 1, (78) 

is split [i is determined by Eqn (J69J). There are also problems related to the topology of the flow diagram. According 



to Fig. \T7\ transitions where quantum number i changes by more than unity are impossible [3J3, ]36j] . However, when 
theorists interpret many experimental data, they state that such transitions occur (see, e.g., review j33| and the 
references therein). 

The general theory of quantum phase transitions does not consider the problem of the relative position of various 
transitions in the phase plane. This theory describes the critical vicinity of one certain transition. Proceeding from 
the assumptions that (a) a transition exists and (b) the factorization £ = £(5x), L v — L^T) occurs in its critical 
region, the resistivity can be described by Eqns (|34]l - l(36|) (see Sections 2.3 and 5.3). Then, the transition point x = x c 
and critical indices can be determined by processing the R(x, T ) curves. This procedure was performed in paper [33| . 
which was discussed in Section 6. The value of the results obtained in [33| becomes clear after the questions formulated 
above and following from flow diagrams are answered. 
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